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An  approximate  numencal  procedure  has  been  developed  for  rapidly  evaluating  the  spectral 
density  integral  that  predicts  the  output  of  a  planar  array  of  many  sensors  excited  by  turbulent 
boundary  layer  pressure  fluctuations.  This  procedure  is  particularly  useful  in  cases  where  the 
transfer  function  factor  of  the  integrand  is  not  a  simple  function  of  the  wave  numbers  in  the  flow 
and  transverse  directions.  The  procedure  exploits  the  facts  that  the  entire  integrand  is  a 
separable  function  of  these  two  wave  numbers  and,  when  the  number  of  sensors  is  large,  the 
array  function  factor  of  the  integrand  is  a  rapidly  varying  function  of  wave  number, 
characterized  by  many  similar  shaped  lobes.  In  addition,  a  model  for  multilayered  media  is 
employed  to  provide  the  transfer  function  for  boundary  conditions  that  closely  correspond  to 
reality.  Results  generated  by  this  procedure  were  compared  to  those  from  an  exact  evaluation  of 
the  integral  which  is  possible  if  the  transfer  function  is  taken  to  be  constant;  there  was  agreement 
to  within  0.2  dB  or  better  over  a  broad  frequency  interval.  Some  results  for  a  realistic  transfer 
function  are  presented,  such  as  the  case  of  an  elastomeric  layer  backed  by  an  elastic  ;,1ate  with 
the  sensors  embedded  at  an  arbitrary  position  within  the  layer. 


INTRODUCTION 

In  a  recent  article,  Ko  and  Nuttall'  discuss  the  evalu¬ 
ation  of  the  spectral  density  integral  that  predicts  the  re¬ 
sponse  of  a  planar  array  of  identical  sensors  when  excited 
by  boundary  layer  pressure  fluctuations  associated  with 
turbulent  fluid  flow  over  the  array.  It  is  assumed  that  the 
turbulent  field  is  spatially  homogeneous  and  temporally 
stationary.  This  integral,  denoted  Q(q})  where  w  is  the 
radian  frequency,  has  the  form 

<2(w)=2ir  r  r  A{k^,ky)S{k^,ky)T{,k^,ky,(o) 

J  —  00  J  —  95 

XP(k^,ky-,a)dk^dky,  (1) 

where  P{k^,ky\o))  is  the  turbulent  boundary  layer  pressure 
spectrum.  The  transfer  function  T{k,^,ky\(o)  describes  the 
wave-number  filtering  action  of  a  layer  of  material,  typi¬ 
cally  an  elastomeric  material,  applied  over  the  sensors  to 
offset  the  turbulent  boundary  layer  from  the  sensors.  The 
transfer  function  expresses  the  magnitude  of  the  ratio  of 
the  pressure  at  any  point  within  that  layer  of  material  to 
the  stress  imposed  by  the  pressure  of  the  turbulent  bound¬ 
ary  layer  acting  on  the  “wet”  face  of  the  layer.  The  hydro¬ 
phone  function  S(k^,ky)  is  a  measure  of  the  spatial  aver¬ 
aging  associated  with  the  size  of  each  sensor.  The  array 
function  is  also  a  measure  of  the  spatial  averaging 

effect  associated  with  more  than  one  sensor  in  both  the  .r 

'Portions  of  Chis  paper  were  presented  at  the  1 20th  meeting  of  the  Acous¬ 
tical  Society  of  America  (J,  Acoust.  Soc.  Am.  Suppl.  1  88.  S6  ( 1990)], 


and  y  directions.  Finally,  k^  and  ky  are  the  wave  numbers 
in  the  x  and  y  directions,  respectively.  The  array  of  sensors 
lies  in  the  x,  y  plane  and  the  turbulent  flow  is  assumed  to 
be  in  the  x  direction.  It  should  be  noted  that  the  expression 
for  Q{(o)  in  Eq.  ( 1 )  appears  to  differ  by  a  factor  of  2ir  from 
that  given  in  a  previous  article  by  Ko  and  Schloemer.^ 
There  is,  in  fact,  no  numerical  discrepancy  between  these 
two  expressions  because  the  factor  of  In  has  been  incorjx)- 
rated  into  the  integrand  in  the  expression  stated  in  Ref.  2. 

It  is  necessary  to  examine  the  functional  form  of  the 
factors  of  the  integrand  in  Eq.  ( 1 )  to  both  appreciate  the 
problem  at  hand  and  to  set  the  stage  for  the  method  that  is 
espoused  to  evaluate  the  integral.  Consider  first  the  turbu¬ 
lent  boundary  layer  pressure  spectrum,  P(k^,ky:a)).  Vari¬ 
ous  analytic  descriptions  of  this  spectrum  have  been  pro¬ 
posed.  A  classic  one,  and  that  which  was  used  by  Ko  and 
Nuttall,'  is  due  to  Corcos.^  Recently,  Ko'*  also  considered 
a  modification  that  predicts  a  much  lower  spectrum  level  at 
very  small  wave  numbers  than  the  original  Corcos  model. 
Because  of  the  desire  to  compare  results  with  those  of  Ref. 
1,  only  the  original  Corcos  model  is  considered  here; 

^(0,w)(a|a2 /Cr) 

(2) 

where  k^—a/u^.  is  the  convective  wave  number  and  is 
the  convective  flow  speed  that  is  taken  to  be  60%  of  the 
free-stream  velocity  U.  The  flow  is  assumed-  to  be  only  in 
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the  X  direction.  The  quantities  a,  and  02  are  dimensionless 
constants  that  are  chosen  to  efTect  a  good  match  of  Eq.  (2) 
with  data.  The  set  of  values  used  in  Ref.  1  is  a,  =0.09  and 
02=0.63.  The  quantity  F(0,6j)  is  called  the  point  power 
spectrum  and  is  described  by 

P{  0,(u )  =  aiplv^/o),  ( 3 ) 

where  pq  is  the  density  of  the  fluid  medium;  is  the 
friction,  or  shear  velocity,  which  is  taken  to  be  3.5%  of  the 
free-stream  velocity;  and  03  is  another  dimensionless  con¬ 
stant  for  which  an  appropriate  value'  is  1.064.  An  impor¬ 
tant  feature  of  the  Corcos  spectrum  as  described  by  Eq. 
(2)  is  that  it  is  a  separable  function  of  the  two  wave  num¬ 
bers  and  ky,  i.e.,  it  can  be  written  as  the  product  of  a 
function  of  kj^  times  a  function  of  ky . 

For  a  rectangular-shaped  sensor  of  dimensions  and 
Ly,  whose  edges  are,  respectively,  oriented  along  the  flow 
and  transverse  directions,  the  hydrophone  function  is^  * 


(%m{k^Ljl)  sin(^3,L/2)\^ 
\  KL,/1  kyLy/2  j  ■ 


(4) 


This  expression  is  simply  the  square  of  the  magnitude  of 
the  normalized  directivity  function  associated  with  the  size 
and  shape  of  the  sensor,  expressed  as  a  function  of  the 
variable  wave  numbers  k^  and  ky  rather  than  as  a  function 
of  the  frequency,  where  the  signal — the  turbulent  pressure 
wave — impinges  upon  the  sensor  at  grazing  incidence.  For 
this  shape  and  orientation  of  the  sensor,  Sikj^.ky)  is  also  a 
separable  function  of  kj,  and  ky.  This  feature  would  not 
exist  for  some  other  shapes  such  as  circular  sensors. 

For  a  rectangular  array  of  equally  spaced,  in  phase, 
identical  sensors,  the  array  function  is’ 


f  sin  ( NJCjd^/2 )  sin  ( NJi/iy/2 )  '\  ’ 
A(k^,ky)  =  sin(My2)  Ny  sin(*//2) )  ’ 

where  and  Ny  are  the  numbers  of  sensors  along  the  x 
and  y  directions,  respectively,  and  and  dy  are  the  center- 
to-center  spacings,  respectively  (Fig.  1 ) .  When  N^  or  Ny  is 
large,  the  array  function,  and  hence  the  integrand  in  Eq. 
( 1 ),  are  rapidly  varying  functions  of  wave  number.  In  such 
cases  conventional  quadrature  techniques  for  performing 
the  integrations  in  Eq.  ( 1 )  can  be  inefiicient  or  impractical. 
This  paper  presents  an  alternative  technique  that  is  both 
fast  and  accurate  for  cases  where  conventional  quadrature 
techniques  become  impractical.  The  method  is  not  appli¬ 
cable  to  small  arrays  {N^  or  A^<5),  however  there  is  no 
need  in  such  cases  because  conventional  methods  may  be 
easily  implemented.  The  array  function  is  the  square  of  the 
magnitude  of  the  normalized  directivity  function  of  an  ar¬ 
ray  of  point  sensors  with  the  same  relative  amplitudes, 
phases,  and  center-to-center  spacings  as  those  of  the  actual 
sensors,  expressed  as  a  function  of  variable  wave  numbers 
k^  and  ky,  again  for  the  case  of  grazing  incidence  excita¬ 
tion.  The  product  of  the  two  functions  A{k^,ky)  and 
S{k^,ky)  is  a  statement  of  the  classic  product  theorem  for 
the  response  function  of  an  array  of  transducers  of  identi¬ 
cal  size,  shape,  and  orientation.  In  the  case  where  beam 
steering  has  been  implemented  to  direct  the  maximum 


FIG.  !.  Plarax  grid  of  rectangular  iV,  along  x  axis  by  Ny  along  y 

axis. 


acoustic  response  of  the  array  to  some  direction  other  than 
broadside,  the  arguments  k^  and  ky  in  the  array  function 
will  accordingly  be  translated  in  wave  number  space  by 
appropriate  amounts  k^x  and  k^y  so  that  the  array  function 
takes  the  form: 


A(kx,ky) 

lsm\Nx{kx-ksx)dx/2]  sin\Ny{ky-ksy)dy/2]  ^ 

~\NxSml(kx-ksx)dx/2]  Ny%in[(ky-ksj,)d/2]j  ' 

(6) 

It  is  noted  that  the  array  function  is  also  a  separable  func¬ 
tion  of  kx  and  ky. 


I.  STATEMENT  OF  THE  PROBLEM 

For  the  special  case  where  7"=  1.0,  i.e.,  there  is  no 
elastomeric  layer  over  the  array  of  sensors,  the  integral  in 
Eq.  (1)  can  be  evaluated  exactly.  Ko  and  Nuttall'  do  this 
both  by  residue  calculus  and  by  a  clever  application  of 
Parseval’s  theorem.  Since  all  three  remaining  factors  in  the 
integrand  of  Eq.  ( 1 )  are  separable  functions  of  kx  and  ky, 
the  double  integral  simplifies  to  the  product  of  two  similar 
integrals,  one  over  the  domain  of  kx  and  the  other  over  ky . 

For  the  case  of  a  realistic  overlayer  between  the  sen¬ 
sors  and  the  turbulent  flow,  however,  the  transfer  function 
T  becomes  complicated.  Ko  and  Schloemer’  discuss  the 
case  of  a  single  homogeneous  elastomer  layer  with  a  rigid 
backing  and  with  the  sensors  embedded  at  an  arbitrary 
position  within  the  layer.  In  this  case,  evaluating  T  in¬ 
volves  solving  a  fourth-order  system  of  linear  algebraic 
equations  whose  coefficients  are  not  simple.  For  a  more 
realistic  condition,  such  as  the  elastomer  layer  on  an  elastic 
base  plate,  T  is  even  more  complicated. 

An  analytical  method  exists  that  permits  rapid  numer¬ 
ical  evaluation  of  the  transfer  function  for  a  harmonic 
plane  wave  incident  upon  a  laterally  infinite  medium  that 
contains  any  number  of  material  layers  bounded  on  either 
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FIG.  2.  Transverse  cross  section  of  multilayered  model  usjd  to  obtain  the 
transfer  function.  Layers  are  laterally  inflnite  and  may  contain  solid  or 
fluid  material.  The  sensor  is  located  at  interface  s. 


side  by  fluids  extending  to  infinity  as  illustrated  by  Fig.  2. 
This  method,  originally  expounded  by  Brekhovskikh,*  is 
based  upon  matching  velocity  and  pressure  variables  at 
each  boundary  of  the  layered  medium.  This  leads  to  an 
analytical  matrix  expression  as  follows: 


where  A'‘  is  a  4x4  matrix  that  relates  velocities  u,,  and 
stresses  at  interface  /  to  the  corresponding  velocities 

and  stresses  at  interface  /+1.  Analytical  expressions  for 
the  components  of  A^*'  are  conveniently  tabulated  in  Ref. 
9.  For  the  purpose  of  computing  a  transfer  function,  one 
needs  a  relationship  that  expresses  the  normal  stress  at 
the  sensor  in  terms  of  the  normal  stress,  ‘  ’,  at  interface 
n+1  where  the  acoustic  plane  wave  of  amplitude  /?q  is 
assumed  to  be  incident.  This  can  be  accomplished  by  ap¬ 
plying  Eq.  (7),  in  turn,  to  each  of  the  intervening  layers 
that  separate  the  sensor  from  the  interface  where  the  exci¬ 
tation  is  incident.  In  this  way,  the  dynamic  variables  at  any 
interface  are  related  to  the  corresponding  dynamic  vari¬ 
ables  at  any  other  interface  by  a  matrix  that  is  the  product 
of  the  A^'*  matrices  for  the  intervening  layers,  the  order  of 
this  product  being  the  same  as  the  order  of  the  intervening 
layers.  Referring  to  Fig.  2,  the  transfer  function  can  be 
expressed  as 


•  (8) 

1  *ZZ  I 

By  imposing  the  appropriate  boundary  conditions  at  inter¬ 
faces  1  and  n -I- 1,  one  obtains  a  set  of  equations  from  which 
T  can  be  evaluated.  Note  that  =  because  the 

fluids  do  not  support  shear  stress.  Also,  the  normal  stress 
at  a  fluid  interface  can  be  related  to  the  normal  velocity;  for 
example,'”  at  interface  n+\ 

(9) 


FIG.  3.  Transfer  function  at  100  Hz  for  a  0,0762-m  (3-in.)-thick  elas- 
tometer  laytr  (material  properties  as  described  in  text)  for  three  different 
base  structures  as  a  function  of  wave  number.  The  sensors  are  located 
0.00635  m  (0.25  in.)  above  the  base. 


where  k  is  the  acoustic  wave  number  in  the  fluid  and  kf  is 
the  wave  number  of  the  imposed  wave  component.  The 
term  containing  accounts  for  the  pressure  radiated 

at  n-(-l,  whereas  Pq  represents  a  spectral  component  of 
flow  noise.  In  this  paper  it  is  assumed  that  fluid  1  is  a 
vacuum  and  therefore  t^'’=0.  Now,  using  the  multilay¬ 
ered  matrix  model  described  above,  one  can  write 


' -f  F33rL'' ^  , 
=  0  =  0414"^  ' ^  +  04211^'’^’^  +  C43rl; 


r‘i'=0=Cj,yi''+”  +  ej2yi''+"-Fe33TL"^‘’.  dOc) 

where  the  components  of  the  P  and  Q  matrices  are  found 
by  taking  the  products  of  the  A*''  matrices  of  the  interven¬ 
ing  layers.  The  transfer  function  can  now  be  obtained  from 
Eq.  (10a),  with  eliminated  by  using  Eq.  (9),  and 

with  and  eliminated  using  Eqs.  (10b)  and 

( 10c).  The  resulting  expression  is  analytic  but  too  complex 
to  permit  an  analytic  evaluation  of  Q{(o).  However,  nu¬ 
merical  evaluation  of  T  at  discrete  values  of  the  wave  num¬ 
bers  is  quite  easy. 

Figure  3  illustrates  the  transfer  function  (decibel 
scale)  as  a  function  of  wave  number  at  a  frequency  of  100 
Hz  for  the  case  of  a  0.0762-m  (3-in.)  elastomer  layer  on 
different  base  structures;  i.e.,  a  perfectly  rigid  base  or  two 
different  thicknesses  of  a  vacuum-backed  steel  plate.  The 
parameters  of  the  elastomer  material  are  chosen  as  follows: 
shear  speed,  16.7  m/s;  dilatational  speed,  2000  m/s;  shear 
loss  factor,  0.33;  dilatational  loss  factor,  0.03;  density  1200 
kg/m^.  These  parameter  values  do  not  represent  a  realistic 
material.  They  were  simply  chosen  to  correspond  to  the 
values  used  in  Ref.  2.  The  point  within  the  layer  where  the 
transmitted  pressure  is  monitored  is  taken  to  be  0.00635  m 
(0.25  in.)  above  the  bottom  face  of  the  elastomer  layer. 
Observe  that  the  two  curves  for  bases  plates  of  finite  thick¬ 
ness  are  significantly  different  from  the  curve  for  the  rigid 
base  at  low  wave  numbers.  The  low  wave-number  spike 
seen  in  all  three  curves  occurs  at  the  acoustic  wave  num¬ 
ber.  The  subsequent  large  excursions  on  the  two  different 
steel  plate  curves  occur  at  their  respective  flexural  wave 
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FIG.  4.  The  four  factors  of  the  integrand  of  the  Q(to)  integral  plotted  as 
a  function  of  wave  number  at  a  frequency  of  100  Hz  for  a  10  by  10  array 
of  0.0508-m  (2-in.)-square  contiguous  sensors.  All  other  parameters  as 
described  in  the  text.  The  transfer  function  curve  is  the  rigid  base  curve 
from  Fig.  3.  The  reference  value  for  the  P  curve  and  the  ASTP  curve  is  1 .0 
(Pa-m2/Hz). 


numbers.  All  the  curves  coalesce  for  large  wave  numbers, 
showing  a  rolloff  that  occurs  around  the  shear  wave  num¬ 
ber  in  the  elastomer  layer.  There  is  a  significant  difference 
at  low  wave  numbers  between  the  curves  for  the  rigid 
boundary  condition  and  the  elastic  plate  base  boundary 
condition,  and  this  will  subsequently  be  seen  to  have  a 
significant  effect  upon  the  value  of  2(<a). 

Suffice  it  to  say  that  in  the  absence  of  a  simple  analytic 
expression  for  T,  the  integral  in  Eq.  ( 1 )  cannot  be  evalu¬ 
ated  to  obtain  a  closed  form  expression  for  Qici>).  Further¬ 
more,  standard  numerical  quadrature  routines  have  been 
found  to  require  copious  amounts  of  CPU  time.  This  arti¬ 
cle  discusses  an  approximate  numerical  procedure  for  eval¬ 
uating  the  integral  in  this  circumstance  that  has  been  found 
to  be  both  fast  and  quite  accurate. 

II.  ANALYSIS 

Figure  4  is  a  plot  of  each  of  the  four  factors  that  com¬ 
prise  the  integrand  of  Eq.  ( 1 )  and  their  product  as  a  func¬ 
tion  of  wave  number  in  the  flow  direction  for  the  fixed 
value  of  ky=0  and  for  the  following  choice  of  parameters: 
t/=  10.288  m/s  (20  kn.),  N^=Ny=l0,  L^=Ly=d^=dy 
=0.0508  m  (2  in.),  ks^=ksy=0,  frequency  =  100  Hz,  and 
the  transfer  function  is  the  curve  for  a  rigid  base  boundary 
condition  in  Fig.  3.  A  similar  plot  vs  ky  would  be  obtained 
for  a  fixed  k^,  except  there  is  no  peak  in  the  pressure 
spectrum  for  any  wave  number  in  the  direction  transverse 
to  the  flow.  Because  of  the  product  nature  of  the  integrand, 
Eq.  ( 1 )  reduces  to  the  product  of  two  similar  integrals,  one 
over  the  domain  of  k^,  and  the  other  over  ky.  It  is  evident 
that  the  controlling  factor  among  the  four  is  the  array 
function  A{k^,ky)  in  the  sense  that  it  is  the  most  rapidly 
varying  function  of  the  integration  variable  k^  (or  ky)-,  i.e., 
over  the  range  of  k^  or  ky  that  defines  one  lobe  of  the  array 
function,  the  other  three  factors  S,  T,  and  P  do  not  vary 
greatly  and  can  be  considered  to  be  essentially  constant. 
This  feature  becomes  more  pronounced  with  increasing 
or  Ny  as  more  and  more  sidelobes  are  present  between 
adjacent  grating  lobes  of  the  A  function.  Furthermore,  the 


> 


shape  of  every  one  of  the  sidelobes  of  the  array  function  , 
appears  to  be  the  same  (note  that  the  discretization  used  in 
generating  this  plot  has  blurred  the  fact  that  there  are  true 
zeros  of  the  array  function).  Consequently,  if  the  integral 
over  any  one  sidelobe  of  just  the  array  function,  when 
renormalized  to  a  maximum  height  of  unity,  was  found  to 
be  a  constant,  i.e.,  independent  of  wave  number,  then  the 
integral  of  the  product  ASTP  over  that  sidelobe  could  be 
approximated  as  that  constant  times  the  values  of  the/4,  5, 

T,  and  P  functions  evaluated  at  the  wave  number  corre- 
sptinding  to  the  peak  of  that  sidelobe.  Integration  over  all 
sidelobes  would  then  entail  a  simple  summation. 

Consider  the  integral  of  one  sidelobe  of  the  normalized 
array  function;  call  this  quantity  normalize  a  side¬ 

lobe  to  unit  height,  note  that  the  factor  sin  (k^^/2)  in  the 
denominator  is  not  a  rapidly  varying  function  of  k^  over 
the  interval  of  that  sidelobe.  Hence,  the  quantity  A^sin 
(k^yi)  is  equated  to  unity,  and  the  integral  in  question 
can  be  written 


r2ir/A-, 


[sin(A^)]2t/f=^, 


(11) 


where  ^ = kjdjl.  The  indicated  limits  define  the  first  side¬ 
lobe;  however,  the  result  is  found  to  be  the  same  for  any 
other  normalized  sidelobe.  The  integral  in  the  ky  direction 
over  any  one  sidelobe  is  formally  the  same  result  except  for 
an  obvious  change  in  subscripts. 

The  mainlobe  and  all  grating  lobes  of  the  A  function 
are  of  different  shape  than  the  sidelobes,  and  these  lobes 
extend  over  a  range  of  wave  numbers  that  is  twice  as  large 
as  that  which  defines  a  sidelobe.  The  integral  over  half  of 
the  mainlobe  or  any  grating  lobe,  denoted  /q,  takes  the 
form 


sin 

Jo 

sin  f 

2 


(12) 


This  integrand  is  already  normalized  to  a  maximum  value 
of  unity.  To  evaluate  this  it  is  convenient  to  note  that 


Af^sin^:  NA 

for  Nj,  odd  and  >  1 
7  /vy2 


(13a) 


-  X  cos(2/— 1)^,  for  even  and  >2. 
1=1 

(13b) 


Hence, 


7“ b: 


(A,-i)/2  (A,-i)/2  /sin[(i-y)2ir/iV,] 

+  11 


+- 


1=1  y=i 
sin[(i-t-7)2w/iV;j] 

77+7) 


(/-/•) 

for  odd  and  >  1 
(14a) 
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FIG.  5.  Comparison  of  the  integrand  of  Eq.  ( I )  for  the  cases  r=  1.0  and 
r  as  in  Fig.  4.  All  other  parameters  as  in  Fig.  4. 


2  ''£  Y  /sinf  (i+2 -  1  )2ir/iVJ 


“l 


+- 


sin[(/— y)2ir/A^jt] 

U-j) 


)■ 


(<+y-i) 

for  N;i  even  and  >2. 

(14b) 


This  result  for  the  integral  over  one-half  of  either  the  main- 
lobe  or  any  grating  lobe  was  compared  with  the  result  in 
Eq.  ( 1 1 )  for  the  integral  over  any  normalized  sidelobe.  The 
result  in  Eq.  ( 14)  is  somewhat  smaller  and  for  large 
approaches  a  value  which  is  approximately  90%  of  Eq. 
( 1 1 ).  In  the  computational  procedure  described  below,  the 
exact  value  of  /g  as  computed  from  Eq.  ( 14)  is  used. 

The  integral  in  Eq.  ( 1 )  is  then  approximated  as  fol¬ 
lows.  Reiterating  that  the  integrand  is  a  separable  function 
of  and  ky,  and  that  the  integrals  over  the  kj,  and  ky 
variables  are  similar  in  form,  consider,  therefore,  only  the 
integral 


A{kAS{kAT{kjc;(o)P{kjc;ci})dkj^. 


(15) 


First,  large  finite  values  are  chosen  for  the  upper  and  lower 
limits  based  upon  knowledge  of  how  the  integrand  de¬ 
creases  with  increasing  ]  k;,  j .  The  controlling  factor  in  this 
regard  is  usually  T{k^,a).  Figure  5  shows  the  integrand 
for  two  different  values  of  T\  the  solid  line  curve  repeats 
the  ASTP  curve  from  Fig.  4,  while  the  dashed  line  curve  is 
the  case  7’=  1.0.  Since  the  transfer  function  is  an  even 
function  of  k^,  limits  of  ±150  m“'  would  appear  to  be 
sufficient  for  integration  of  the  former  curve,  while  the 
limits  should  be  increased  to  about  ±300  m~'  for  the 
latter.  Since  the  integrand  is  positive,  a  lower  bound  on  the 
exact  value  will  be  realized  by  restricting  the  limits.  As¬ 
suming  the  functions  S,  T,  and  F  to  be  essentially  constant 
over  the  interval  of  either  one  sidelobe  or  one-half  of  either 
the  mainlobe  or  a  grating  lobe,  QA<^)  **  then  approxi¬ 
mated  as 


Cx(«)=  X-S'(*x,/)^(*x.,;«)^(*x,/;«)  r  A{k^)dk^, 

/  J  lobe, 


(16) 


where  the  indicated  summation  extends  over  all  the  lobes 
within  the  chosen  limits  of  integration.  The  values  k^^  de¬ 
note  the  maxima  of  the  A  function  in  the  interval  of  inte¬ 
gration  and  are  determined  by  finding  the  roots  of  the 
transcendental  equation  that  results  from  setting  the  deriv¬ 
ative  of  the  A  function  with  respect  to  k^  equal  to  zero: 


(17) 


The  remaining  integral  in  Eq.  (16)  is  then  either  approx¬ 
imated  as  Tr/JVjjd^  times  the  value  of  the  array  function  at 
k^j  for  any  sidelobe,  or  the  value  predicted  by  Eq.  ( 14)  for 
either  half  of  the  mainlobe  or  any  grating  lobe.  Hence, 


'NM 


(18) 


where  a,=  l  if  /  denotes  a  sidelobe,  but  a,=A^^g/7r  if  i 
denotes  either  half  of  the  mainlobe,  or  either  half  of  any 
grating  lobe;  again,  the  latter  values  of  a,  vary  with  and 
are  numbers  slightly  greater  than  0.9. 

Since  the  integral  over  ky  is  formally  the  same,  the  final 
result  for  Q(ei>)  is 


277^  /  \ 

g(a))=  -  -  .  (  ^a^(k^i)S(k^j)T(k^i;cj)P(k^j;o}) 


a  jA  ( kyj)S  ( kyi)  T  ( kyj  ;(o )  P{  kyj  ;co ) 


(19) 


III.  RESULTS  AND  DISCUSSION 

A  FORTRAN  language  computer  program  was  devel¬ 
oped  to  generate  numerical  results  from  Eq.  (19).  An 
IMSL  routine  was  used  to  find  the  roots  of  Eq.  ( 17),  i.e., 
the  k^(  values,  in  the  interval  between  the  peak  of  the 
mainlobe  and  the  peak  of  the  first  grating  lobe.  The  re¬ 
mainder  of  the  kji  i  values  within  the  chosen  limits  of  inte¬ 
gration  are  obtained  simply  by  noting  the  periodic  nature 
of  the  array  function.  These  values,  shifted  by  k^x  if  appro¬ 
priate,  are  then  stored.  The  same  procedure  produces  the 
set  kyj.  Note  that  under  this  counting  system  the  entire 
mainlobe  and  every  grating  lobe  is  accounted  for  by  two 
successive  values  of  index  i  or  j  with  identical  values  of  k^  , 
or  kyj.  The  indicated  summations  in  Eq.  (19)  are  then 
performed  in  straightforward  fashion  to  produce  Q(a)). 

The  analysis  and  computer  program  were  checked  by 
calculating  Q{a>)  for  the  set  of  parameter  values  used  by 
Ko  and  Nuttall,'  viz.,  a  10  by  10  array  of  contiguous 
0.0508-m  (2-in.) -square  sensors  with  no  elastomer  over¬ 
layer  ( 7’=  1.0),  a  free-stream  velocity  of  10.288  m/s  and  a 
frequency  range  of  100  Hz  to  2  kHz.  The  results  obtained 
matched  their  exact  values  within  0.2  dB  over  the  whole 
frequency  range.  In  all  cases  the  values  obtained  from  the 
approximate  expression  are  smaller  than  the  exact  values. 
The  approximate  expression  should  become  more  accurate 
as  the  frequency  of  the  sidelobes  of  the  array  function  A 
increases  (i.e.,  as  either  parameter  Nx  and  Ny  increases) 
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Fr«qu9ncy  (Hz) 

FIG.  6.  Spectral  density  level  for  various  situations:  (a)  one  point  sensor 
and  unity  transfer  function;  (b)  a  10  by  10  array  of  contiguous  O.0SO8-m 
(2'in.)-square  sensors  and  unity  transfer  function;  (c)  case  (b)  with  a 
0.0762-m  ( 3-in. )  overlayer  of  elastomer  with  a  rigid  base,  sensors  0.00635 
m  (0.25  in.)  above  base;  (d)  same  as  (c)  except  base  is  a  vacuum  backed 
0.0254-m  ( l-in.)-thick  steel  plate.  All  other  parameters  as  defined  in 
the  text. 

because  the  assumption  that  the  S,  T,  and  P  functions  are 
constant  over  a  sidelobe  of  the  A  function  then  becomes 
more  valid.  That  this  is  the  case  can  be  demonstrated  by 
subdividing  the  array  into  a  larger  number  of  smaller  con¬ 
tiguous  sensors  that  cover  the  20-in.-square  surface.  The 
exact  result  for  Q(<u)  should  be  the  same.  Hence,  the  ap¬ 
proximate  expression  was  used  to  calculate  Qio))  for  a  20 
by  20  array  of  0.0254-m-square  sensors  and  40  by  40  array 
of  0.0127-m-square  sensors.  The  computation  time  was  not 
perceptibly  increased  in  spite  of  the  doubling  and  quadru¬ 
pling  of  the  number  of  terms  to  be  summed.  As  expected, 
the  accuracy  of  approximation  improved.  The  Q(.(o)  values 
increased  slightly  to  within  0. 1  dB  of  the  exact  values. 

To  illustrate  the  effect  of  including  a  nonconstant 
transfer  function  in  the  computations,  the  results  shown  in 
Fig.  6  were  generated.  For  wave  number  tilts  ksx  and  k^y  of 
zero,  the  spectral  density  level  is  plotted  as  a  function  of 
frequency  for  the  cases  of  (a)  one  point  sensor  with  no 
elastomer  overlayer  so  that  r=  1.0;  (b)  a  10  by  10  array  of 
contiguous  0.0508-m-square  sensors,  again  with  T  =  1 
(these  first  two  cases  are  data  from  Ref.  1);  (c)  case  (b) 
with  a  0.0762  m  thick  overlayer  of  elastomer  (material 
properties  defined  earlier)  with  a  rigid  base  where  the  sen¬ 
sors  are  0.00635  m  above  the  rigid  base;  (d)  same  as  (c) 
except  the  rigid  base  is  replaced  by  a  vacuum  backed 
0.0254-m-thick  steel  plate.  The  transfer  functions  for  cases 
(c)  and  (d)  are  plotted  in  Fig.  3  for  a  frequency  of  100  Hz. 
There  is  evidently  not  much  difference  in  Q((x))  for  cases 
(b)  and  (c).  As  is  apparent  from  Fig  3,  the  transfer  func¬ 
tion  at  100  Hz  for  the  rigid-backed-elastomer-layer  situa¬ 
tion  is  essentially  unity  up  to  a  wave  number  of  20  m  “  * 
which,  noting  the  plots  of  the  integrand  in  Fig.  5,  is  large 
enough  to  include  most  of  the  wave  number  region  that 
contributes  to  Q(.ci}).  Consequently,  the  integral  is  not 
much  affected  by  the  high  wave  number  rolloff  of  T.  Con¬ 
versely,  the  smaller  transfer  function  associated  with  the 
more  realistic  case  of  the  vacuum  backed  steel  base  plate 
results  in  a  noticeably  lower  value  of  Q(ci>),  particularly  at 
the  low  end  of  the  frequency  range.  The  bumps  in  curve 


(d)  at  440  and  1600  Hz  are  consequences  of  certain  peaks 
in  the  transfer  function,  such  as  the  one  associated  with  the 
shear  wavelength,  becoming  coincident  with  the  peak  of  a 
high  sidelobe  of  the  array  function.  Note  that  in  this  special 
case  of  contiguous  sensors,  the  nulls  of  the  5  function  co¬ 
incide  with  the  peaks  of  the  grating  lobes  of  the  A  function, 
and  hence  there  are  no  grating  lobes.  Were  that  not  the 
case,  there  would  probably  also  be  a  large  bump  in  the 
Q((i>)  curve  at  frequencies  where  the  peak  of  the  P  curve 
coincided  with  a  grating  lobe. 

iV.  CONCLUSION 

An  approximate  numerical  technique  has  been  devel¬ 
oped  for  evaluating  the  spectral  density  integral  over  all 
wave  numbers  to  obtain  the  output  of  a  large  planar  array 
of  sensors  being  excited  by  turbulent  boundary  layer  pres¬ 
sure  fluctuations.  This  technique  is  fast  [for  comparison, 
Q((i))  for  the  above  case  where  7"=  1.0  was  also  evaluated 
using  the  IMSL  numerical  quadrature  routine  QDAGS 
and  the  elapsed  CPU  time  was  found  to  be  almost  twice  as 
much  as  for  this  approximate  technique]  and  is  quite  ac¬ 
curate  as  determined  by  comparing  the  results  with  those 
for  cases  where  the  integrals  in  question  can  be  evaluated 
exactly.  Most  importantly,  it  can  be  used  when  it  is  not 
possible  to  evaluate  the  integrals  exactly,  such  as  when  a 
simple  analytic  expression  for  the  transfer  function  is  not 
available.  All  that  is  required  for  this  technique  is  knowl¬ 
edge  of  the  value  of  the  transfer  function  at  a  discrete  set  of 
wave  number  values.  It  has  also  been  demonstrated  how 
numerical  values  for  the  transfer  function  may  be  easily 
obtained  for  multilayered  array  models  that  more  closely 
resemble  reality  than  for  models  that  assume  rigid  sub¬ 
strates.  The  technique  for  computing  Q(a))  for  exploits  the 
fact  that  the  array  function  is  the  most  rapidly  varying 
function  of  wave  number  of  all  the  factors  of  the  integrand, 
and  assumes  that  the  other  factors  may  be  considered  to  be 
constant  over  the  wave  number  interval  corresponding  to 
one  lobe  of  the  array  function.  Consequently,  the  accuracy 
of  the  technique  increases  as  the  number  of  sensors  in¬ 
creases,  since  that  number  is  a  measure  of  the  number  of 
lobes  between  adjacent  grating  lobes  of  the  array  function. 
The  technique  is  restricted  to  situations  where  all  the  fac¬ 
tors  of  the  integrand  are  separable  functions  of  the  wave 
numbers  in  the  flow  direction  and  transverse  to  that  direc¬ 
tion.  Hence  a  regular  grid  of  circular  shaped  sensors,  or  an 
irregularly  spaced  array  of  rectangular  sensors,  could  not 
be  evaluated  by  this  technique.  The  technique  was  illus¬ 
trated  for  the  case  of  contiguous  sensors  but  is.  in  no  way, 
limited  to  this  case.  Gaps  of  any  size  may  be  considered 
between  sensors  as  long  as  the  gap  sizes  are  uniform  along 
each  lateral  direction. 
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